STABILITY OF SOLUTIONS TO NONLINEAR 
DIFFUSION EQUATIONS 



TEEMU LUKKARI 

Abstract. We prove stability results for nonlinear diffusion equa- 
tions of the porous medium and fast diffusion types with respect 
to the nonlinearity power m: solutions with fixed data converge in 
a suitable sense to the solution of the limit problem with the same 
data as m varies. Our arguments are elementary and based on a 
general principle. We use neither regularity theory nor nonlinear 
semigroups, and our approach applies to e.g. Dirichlet problems in 
bounded domains and Cauchy problems on the whole space. 



1. Introduction 
We study the stability of positive solutions to the parabolic equation 
(1.1) d t u - Au m = 

with respect to perturbations in the nonlinearity power m. The main 
issue we address is whether solutions with fixed boundary and initial 
data converge in some sense to the solution of the limit problem as 
m varies. This kind of stability questions are not only a matter of 
merely mathematical interest; in applications, parameters like m are 
often known only approximately, for instance from experiments. Thus 
it is natural to ask whether solutions are sensitive to small variations 
in such parameters or not. 

The equation (11. ip is an important prototype of a nonlinear diffusion 
equation. For m > 1, this is called the porous medium equation (PME), 
and m < 1, the fast diffusion equation (FDE). The PME is degenerate, 
the diffusion being slow when u is small. The FDE is singular, and 
the opposite happens: the diffusion is fast when u is small. We do not 
exclude the case m — 1, when we have the ordinary heat equation. 
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However, we do restrict our attention to the supercritical range: 

m > m c , where m c = in — 2) + /n. 

For the basic theory of the porous medium and fast diffusion equations, 
we refer to the monographs |6j [26], |27] and the references therein. 

It turns out that stability with respect to m for (II. ip can be es- 
tablished in a relatively elementary manner. The principal reason for 
this is that weak solutions to (11. ip are defined in terms of the function 
u m instead of u. This means that u m and its gradient are always L 2 
functions, even if m varies. The situation is markedly different for the 
]9-parabolic equation, 

(1.2) d t u - divfl Vu\ p ~ 2 Vu) = 

where changing the exponent p also changes the space to which weak 
solutions belong. Dealing with this in the case of (11.2]) is quite delicate, 
and one needs the fact that the gradient satisfies a reverse Holder in- 
equality. The equation (11.11) has the advantage that we may work in 
a fixed space. For this reason, our results do not rely on the more so- 
phisticated tools of regularity theory. More spesifically, we do not need 
Holder continuity, Harnack's inequality, or reverse Holder inequalities 
for the gradient. 

The starting point of our argument is a compactness property of 
weak solutions: locally uniformly bounded sequences of weak solutions 
contain pointwise almost everywhere convergent subsequences. With 
the compactness result in hand, a stability result with convergence in 
an LP space for a particular problem follows by verifying the local uni- 
form boundedness and that the correct initial or boundary values are 
attained. Only some fairly simple estimates are needed for the sec- 
ond step. We carry out this step in detail in two cases: for Dirichlet 
problems with nonzero boundary values on bounded domains, and for 
Cauchy problems on the whole space, with initial data a positive mea- 
sure of finite mass. In the latter case, our stability result applies to the 
celebrated Barenblatt solutions. 

The previous result closest to ours is that of Benilan and Crandall [2]. 
They prove the stability of mild solutions to 

(1.3) d t u - A<p(u) = 

with respect to if using the theory of nonlinear semigroups. Here ip can 
be a maximal monotone graph; in the case of a power function, the 
assumption used in [2j reduces to m > m c . Our approach is different 
from that of [2] , as we do not employ the machinery of nonlinear semi- 
groups. Explicit estimates covering the situation of [2J are established 
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in pi]. See also |21[ 122] for the onedimensional situation. Stability re- 
sults are also known for, e.g., equations similar to the p-Laplacian |17| . 
the p-parabolic equation (II. 2p |15| . obstacle problems [16], and eigen- 
value problems [18]. See also [3j 123] for the limit when m — > oo, and 
[101 K21 OS] for m -> 0. 

Generalizing our argument to other boundary and initial conditions 
is straightforward. For instance, Neumann boundary conditions and 
Cauchy problems with growing initial data can be handled similarly. 
Generalization to nonlinearities other than powers, as in (II. 3p . should 
also be possible, albeit less straightforward. 

Our results include the case when the limiting problem is the heat 
equation. In this situation, we do not need the restrictions m > 1 
or m < 1 on the approximating problems; both degenerate (m > 1) 
and singular (m < 1) problems are allowed. However, the restriction 
m > m c seems essential, as a number of the tools we use are known 
to fail when < m < m c . For instance, local weak solutions might 
no longer be locally bounded, and the L l -L°° smoothing effect for the 
Cauchy problem fails. 

We also provide an alternative argument for the stability of Dirichlet 
problems in bounded domains in the case m > 1. Compared to the 
previous argument, the advantages of this approach are the fact that 
local boundedness is not needed, and that one can in addition estimate 
the difference of two solutions in terms of the difference of the respective 
nonlinearity powers. The disadvantages are the restriction m > 1 and 
the fact that the proof seems less amenable to generalizations, since we 
employ strong monotonicity. See [2~Tj |2~2"] for similar estimates in the 
onedimensional case. 

The paper is organized as follows. In Section [2J we recall the neces- 
sary background material, in particular the definition of weak solutions. 
Section [3] contains the proof of the compactness theorem for locally 
bounded sequences of weak solutions. The actual stability results are 
then established in Sections H] and [51 for Dirichlet problems in the for- 
mer and Cauchy problems in the latter. We finish by presenting the 
alternative proof for Dirichlet problems in Section [6j 

2. Weak solutions 

Let Q be an open subset of IR n , and let < t\ < t 2 < T. We use 
the notation Qt = x (0,T) and — U x (ti,^), where U C 

is open. The parabolic boundary d p U tlt t 2 of a space-time cylinder U tl ,t 2 
consists of the initial and lateral boundaries, i.e. 



d p U tlM = {Ux{t x })\J{dUx[t x M)- 
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The notation U tl .t 2 <s Qt means that the closure C7 tl; t 2 is compact and 

u tut2 c q t - 

We use if x (f2) to denote the usual Sobolev space, the space of func- 
tions u in L 2 (Q) such that the weak gradient exists and also belongs 
to L 2 (Q). The norm of H\Q) is 

= ||w||i2(n) + ||Vm||_l2(q). 

The Sobolev space with zero boundary values, denoted by Hq(Q), is 
the completion of C^°(Q) with respect to the norm of H 1 (Q). The dual 
of ffi(fi) is denoted by H-\Q). 

The parabolic Sobolev space L 2 (0,T; H 1 ^)) consists of measurable 
functions u : Qt [—00,00] such that x 1— > u(x,t) belongs to H 1 (Q) 
for almost all t G (0, T), and 

\u\ 2 + I Vw| 2 dx dt < 00. 

The definition of L 2 (0, T; Hq (Q)) is identical, apart from the require- 
ment that x u(x,t) belongs to Hq(Q). We say that u belongs to 
Ll c (0,T;Hl oc (Q)) Hue L 2 {t u t 2 ; H\U)) for all U tlM m Q T - 

We use the following Sobolev inequality. See Proposition 3.1, p. 
7] for the proof. 

Lemma 2.1. Let u be a function in L 2 (0,T; Hq(Q)) . Then we have 

(2.1) f \u\ 2K dxdt<C [ \\7u\ 2 dxdt fesssup f u l+l/m dx] ' , 
Jn T Jn T V o<t<T Jn J 

where 

, , 1 1 
2.2) «=! + _ + . 

n mn 

Solutions are defined in the weak sense in the parabolic Sobolev 
space. 

Definition 2.2. Assume that m > m c . A nonnegative function u : 
Qt — t- M is a local weak solution of the equation 

f)u 

(2.3) — - Au m = 
in Q T , if u m eL 2 oc (0,T;HUQ)) and 

(2.4) f -u^f + Vu m -V<pdxdt = 

for all smooth test functions tp compactly supported in Qt- 
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We will always assume that m > m c , and consider only nonnegative 
solutions. We refer to the monographs jSJ ESJ EZ] for the basic theory 
related to this type of equations, and numerous further references. In 
particular, weak solutions have a locally Holder continuous representa- 
tive, see [5j[8] or Chapter 7 of [27]; however, we do not use this fact. 

An important example of a local weak solution in the sense of Defi- 
nition [2]2] is provided by the celebrated Barenblatt solution [TJ [28]. For 
m > 1, it is given by 




A(m-l) |x| 

(2.5) B m (x,t) = {° V" 2mn * 2A/n 



2 \l/(m-l) 

i 

t < 



, t>0, 



where 



A 



n 



n{m — 1) + 2 

For m c < m < 1, it is convenient to write the formula as 

( t -*(c+k Nir 1/{1 - m) t>0 

(2.6) B m (x,t)=r \ U + K t^) > t>u > 

0, t < 0, 



where 



A(l -m) 



2/11/1 

Note that A and k are strictly positive, since here m c < m < 1. 

Observe that for m > 1, £> m has compact support in space at each 
time instant, while for m < 1, £> m has a powerlike tail. In both cases, 
B m is a local weak solution in IR n x (0, oo). See Section [5] for the trace 
of B m at the initial time t = 0. 

The definition of weak solutions and supersolutions does not include 
a time derivative of u. However, we would like to use test functions 
depending on u, and thus the time derivative ^ inevitably appears. To 
deal with this defect, a mollification procedure in the time direction, for 
instance Steklov averages or convolution with the standard mollifier, is 
usually employed. The mollification 

i r* 

(2.7) u*(x,t) = - e {s - t)/a u(x,s)ds 

° Jo 

is convenient. The aim is to obtain estimates independent of the time 
derivative of u*, and then pass to the limit a — > 0. 

The basic properties of the mollification (12. 7p are given in the fol- 
lowing lemma, see 1201 . 



6 TEEMU LUKKARI 

Lemma 2.3. (1) If u E L p (Vt T ), then 

||«*||i>(n T ) < ||w|Up(n r ), 

du* u — u* 



2.8 



8t ~ a ' 

and u* — >■ w m L p (fi*r) as cr — ^ 0. 

(2) //Vm G L p (fi T ), t/ien V(u*) = (Vu)*, 

||Vw*||iP(n T ) < || V-u||LP(n T ), 

and Vu* — >■ Vu m L p (Vt T ) as a — > 0. 

(3) Ifuk^-u in L p (Qt), then also 

. Out du* 
u k -)■ u* and -)> 



at at 



zn L p (fi T ). 

(4) //Vu* -> Vu zn L P (Q T ), then Vu* k ->■ Vu* zn L P (Q T ). 

(5) Similar results hold for weak convergence in L p (Qt)- 

(6) If f E C(Pr) } then 

<p*(x, t) + e~ t/a <p(x, 0) -> t) 

uniformly in VLt as a — > 0. 

We use the following estimate for the local version of our stability 
result. See [HI Lemma 2.15] or [19, Lemma 2.9] for the proof. 

Lemma 2.4. Let u be a weak solution such that < u < M < oo, 
where M > 1, and let r\ be any nonnegative function in Cg°(fi). Taen 

/ n 2 |Vw m | 2 dxdt < 2M m+1 / n 2 dx + 16M 2m / |Vr/| 2 dxdt. 
Jq t Jn Jq t 

We use the following elementary lemma to pass pointwise conver- 
gences between various powers. 

Lemma 2.5. Let be a sequence of positive functions on a measur- 
able set E with finite measure such that 

fi—tf in L 1 (E)and pointwise almost everywhere. 

Assume that cti^r a as i — > oo. Then 

ft -> r 

pointwise almost everywhere. 
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Proof. By Egorov's Theorem, for any e > there is a set F £ such that 
\F £ \ < e and /j -> / uniformly in E\F e . Pick any point x <E E \ F e 
such that 

< <5 < /(x) < M < oo. 

Then 

< 5/ 2 < /i(^) < 2M 
for all sufficiently large % by uniform convergence. An application of the 
mean value theorem to the function a^f gives 

\f t ( X r - f(xr\ <\fi{xr - Mxr\ + \f t ( X r - /(zri 

<c(5, M)\a t -a\ + \f t (x) a -f(xy\. 

By the convergence assumptions, it follows that fi(x) ai — > f(x) a as 
i — > oo. 

Since 5, M, and e are arbitrary, the above implies that fi(x) ai — > 
f(x) a for almost all x in the set {0 < f(x) < oo}. For points where 
f(x) = 0, it is easy to check that fi(x) at — > 0. Hence we have the 
desired convergence almost everywhere in the set {f(x) < oo}, which 
is sufficient since / is integrable. □ 

3. Stability of local weak solutions 

In this section, we establish a local version of stability for a bounded 
family of local weak solutions. This is beneficial since we may then ap- 
ply the same result to both initial-boundary value problems in bounded 
domains and Cauchy problems on the whole space. The crucial point 
in our stability results is extracting pointwise convergent subsequences 
out of a sequence of solutions, in other words, a compactness property 
of solutions. 

Recall that we use the notation U d to mean that the closure U 
is compact and contained in Q. In view of applying this result to the 
Cauchy problem on IR n , we allow the cases f2 = R n and T = oo. We 
denote 

m + = sup m 8 and m~ = infmj, 

i i 

with similar notations for other exponents. 

The main result of this section is the following theorem. 

Theorem 3.1. Let m^, i — 1, 2, 3, . . . , be exponents such that 

(3.1) mj-^m as i — > oo for some m > m c = (n — 2) + /n 

Let u i; i = 1,2,3,..., be positive local weak solutions to 

d t u t - An™* = 
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in Qt- Assume that we have the bound 

( 3 - 2 ) IKIU °(fft 1 ,t a ) ^ M < 00 

in all cylinders U tl .t 2 such that U d Q and < t% < £2 < 7\ 

TTien £/iere is a function u such that Ui — )■ w and w mi — )■ w m pointwise 
almost everywhere in Qt, for a subsequence still indexed by i. Further, 
u is a local weak solution to 

d t u - Au m = 0. 

We split the proof of Theorem 13.11 into several lemmas. In view of 
the convergence assumption (I3.ip . we are free to assume that 

m c < m~~ and m + < 00. 

Let us define the auxiliary exponents 

m\ = max{mj, 1} and m\ = min{mj, 1}. 

Then mjj > ml, and one of these exponents always equals mi and the 
other equals one. We use a similar notation for the limit exponent m. 

The proof of Theorem 13.11 consists of three main steps: first we apply 
a compactness result j25] to certain auxiliary functions to find the limit 

function u. Then we show that u { ' converges to u m in measure. This is 
the most involved part of the proof, and for a key estimate we apply a 
test function due to Oleinik. In the last step, we establish the pointwise 
convergences and w™ -4 — > u m by applying Lemma [2.51 and the 

fact that u is a local weak solution by applying Lemma [2.41 

The next lemma will be used in proving the convergence in measure. 

Lemma 3.2. Let < s < t < M, where M > 1, be such that 



Then 



£ m * - s mS > A > 0. 



£ mb - s mb > —M m "- m " min{A, A mb/m *} 



Proof. Consider first the case s = 0. Then 

^m? g t m ^ ^> A m V m " 

since £ m " > A. Assume then that s > 0. By the mean value theorem, 
we have 



STABILITY OF SOLUTIONS TO NONLINEAR DIFFUSION EQUATIONS 9 
for some f G (s m \t m *). Since ^ - 1 < and t m * < M m \ we have 



which completes the proof. □ 

The following lemma provides the key estimate for showing that the 
original sequence also converges to the limit found by applying the 
compactness result. 

Lemma 3.3. Let U andui be as in Theorem \3.1\ and let < t\ < t 2 < 
T. 

Fix a number < 5 < I, and suppose that u^s is the unique function 
which satisfies 

d t u i: s - Am™ 1 =0 in U tljt2 , 

Ui t s{x, ti) = ui(x, tx) + 5, x G U. 



Then 



where 



(u i>3 - Ui)(u% - O dxdt < c(5 + 5 m ~), 

U tl ,t 2 



c = 2 m++1 (M m+ + M + l)\U tu 



,t 2 \ 

and M is the number appearing in (13. 2p . 

Before proceeding with the proof let us note that the technical reason 
for introducing the exponents mjj and m- is that we may write 

{ui, s - Ui)K* - u i ) = ( u iJ - U i )( U iJ - u i ) 
when applying this lemma. 

Proof. The proof is an application of a test function due to Olemik. 
The function u™j — n™ - * — 5 mi has zero boundary values in Sobolev's 
sense, so the same is true for the function 



T](x,t) 



f t t2 u™i-ur-5^ds, tl <t<t 2 , 
0, t > t 2 . 
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We use r\ as a test function in the equations satisfied by u^s and U{ 
and substract the results. This gives 



U H ,t 



+ 



v(«3 - < 



<2 



V(m^ — m^) ds dx dt 



--S 




■it 



[7 ./ti 



i,5 



M„ mi ds dx. 



We move the term with 5 to the right hand side, and integrate with 
respect to t in the elliptic term. We get 



/ (u h5 - tti)(«3 - uT) dxdt + \ [ [ f 2 V 



n 2 



«r 4 ) ds 



dx 



= 5 m > / - Mi ) dx dt + 5 / (u™Z - uf*) dx dt 

Ju tl ,t 2 Ju tl ,t 2 

Now, both of the terms on the left are positive, so we have the freedom 
to take absolute values of the right hand side. The claim then follows 
by discarding the second term on the left hand side, and estimating the 
integrals on the last line by using (13. 2p and the fact that Ui t s < M + 1, 
by the comparison principle. Indeed, we have 



(uls - Ui) dxdt 



<2(M + l)\U tut2 \ 



and 



«5 -up) dxdt 



<2 m+ (M m+ + l)\U tut2 \. 



□ 



The following lemma is the key step in the proof of Theorem 13.11 

Lemma 3.4. There exists a subsequence, still indexed by i, and a func- 
tion u such that 

u i — >• u as i — > oo 
in measure and pointwise almost everywhere in U tlt t 2 - 

Proof. Recall that U is an open set such that U <s fi, and < t% < ti < 
T. We use the auxiliary functions Ui t $ defined in Lemma 13.31 for any 
< 5 < 1, Ui t s is the weak solution in Ut x ,t% with boundary and initial 
values up+8 mi and u(-, t L ) + 5, respectively. Note that 5 < u it s < M+5 
by the comparison principle. 
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Fix then a regular open set V (e U, and t\ < S\ < S2 < £2- For 
indices such that rrti > 1, we have 

iv^i =v«^) 1/m * = k 1/mi_1 v<^| 
<8 llmi - x \vu™i\ < <y 1/m+-1 |Vit$|, 

and for indices such that m ; < 1 we have 

|v^| =|v(«3) 1 /^| < i^-Vu^i 

<(M + 5) 1 / m t- 1 |Vu^| < {M + l) 1 ' m ~- 1 \Vu^\. 

Thus Lemma EH implies that (Vu^) is bounded in L 2 (V^ 1)S2 ). The fact 
that the sequence (d t Ui t s) is bounded in L 2 (si, S2] H^iV)) follows easily 
from the equation satisfied by Ui, S and Lemma [2.41 An application of 
[25| Corollary 4] shows that the sequence {u^s) is compact in L 2 (V SltS2 ). 
Thus for any fixed 5 > 0, we may extract a subsequence which converges 
in L 2 (V SlyS2 ) and pointwise almost everywhere. 

The next step is a repeated application of the compactness estab- 
lished in the previous step. Let 8j = 1/j, j = 1, 2, 3, . . . , and for j = 1 
pick indices i such that 

for some function u$ 1 , the convergence being in L 2 (V Sl-S2 ) and pointwise 
almost everywhere. We proceed by picking a further subsequence so 
that also 

Ui,S 2 -> u &2 

in L 2 (V sltS2 ) and pointwise almost everywhere. We continue this way, 
and get a twodimensional table of indices where the values of % 

in the diagonal positions have the property that 

Ui,Sj ~ > us 3 as i — > 00 for all j, in L 2 (l^ ljS2 ) and pointwise a.e. 

By the comparison principle, we have 

Ui,8 3 > Ui,8 k for k > j, 

and letting i tend to infinity we get 

w<5 3 > us k , whenever > j. 

Thus we may define the function u in the claim of the current lemma 
as the pointwise limit 

u(x, t) = lim UgAx, t). 

j-s>oo 3 

We have now found the limit function u, and the proof will be com- 

pleted by showing that u i ' converges to u m in measure, first in V Sl:S2 
and then in U tl j 2 by an exhaustion argument. The first convergence is 
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a consequence of Lemmas 13.21 and 13.31 Obviously — > M m " pointwise 
almost everywhere as j — > oo. Further, Lemma 12.51 implies that also 

Uj x U T a.e. as i — >• oo. 

1,0 j Oj 

To show the convergence in measure, let A > and e > 0. We have 
\{\uf-u m> \>\}\<\{\uf-u±>\/3}\ 



+\{\<X -<l > A/3} I + |{|< > A/3} | 

Since — w™* > cmin{A, A m » //m »} by Lemma |3.2[ we have 

|{|«r ! -^l>A/3}|= / m) mt ldxdt 

J\\u. l -u. } |>A/3> 



3 



< r — 5— / (uTx — ) (uTx — uT* ) dx dt 

Amin{A,A m */ m *} 3 
where the last inequality follows from Lemma 13.31 and the fact that 

{Ui, 5 - Ui){u^ - <*) = { Ui J - Ui *)K« - u, •). 
Thus for each fixed A > 0, we may choose Sj small enough, so that 

tt tt 
{{{UiJ.-v.i '|>A/3}|<e. 

Here it is crucial that this choice can be made independent of i. 

From the pointwise convergence established above, it follows that 
U T ~~ >" i n measure. Thus we may choose 5j small enough, so that 

|{|<"- M m "|>A/3}|<,. 



Finally, we know that u?$. — > Ug 1 . pointwise almost everywhere and 
hence also in measure. Thus the above estimates imply that 



limsup \{\uf - u mS \ > A}| < lim \{\u™} - u^}\ > A/3}] + 2e = 2e. 
Since e was arbitrary, we get 

lim \{\uf -u mi \ > A}| = 
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for all A > 0, as desired. We have proved that for cylinders V SlyS2 such 
that Vg ljS2 d U tl j 2 <s Qt, w e may extract a subsequence such that 

— y u m measure m V SljS2 as % — > oo. 

Passing to a further subsequence, we have convergence pointwise almost 
everywhere. 

To find a subsequence converging in the whole of U tl j 2 , we use a 
diagonalization argument. Exhaust U by regular open sets V , and 
choose nested time intervals (s\, s^) <s (sf , s|) <e . . ., k = 1, 2, . . ., such 
that 

jt 

First, pick a subsequence such that (u i *) converges in measure in a i 

to u m> . The procedure continues inductively, by the selection of a fur- 
ther subsequence that converges in measure in V £\ k+1 to the function 

s 1 ,s 2 

u m . Taking the fcth index in the subsequence selected in the fcth step 
yields a subsequence convergent in measure in U tljt2 . □ 

The next step in the proof of Theorem 13. H is to prove the convergence 
of the other powers of u by using the convergence established in the 
previuos lemma. 

Lemma 3.5. Let Ui, i = 1,2,3,..., be such that < Ui < M < oo, and 

i — > u m in measure in U tlt t 2 - Then Ui — > u and vl™" 1 — > u m almost 
everywhere in U tl j 2 as i oo. 

Proof. The bound Ui < M and convergence in measure imply that 

u™ 1 * — > u m% in L g (U tl j 2 ) f° r an Y finite q. Thus we obtain the desired 
convergences by two applications of Lemma 12. 5\ first passing from the 

convergence u { i — > u m to the convergence Ui — > u, and then from 
Ui -)■ u to n™ 1 -)> u m . □ 

With all of the preceding lemmas available, the proof of Theorem 13. II 
is now a relatively simple matter. 

Proof of Theorem IJ.il The pointwise convergences follow from Lemma 
13.51 by an exhaustion argument, as in the proof of Lemma [3.41 To show 
that u is a local weak solution, fix a test function ip e C^°(Qt)- We 
use the bound Ui < M and Lemma 12.41 to get subsequences of (ui) 
and (Vti^), weakly convergent in L 2 (supp<y?). Due to pointwise con- 
vergences, we see that the weak limits must be u and Vu m , respectively. 



14 TEEMU LUKKARI 



From the weak convergences, it follows that 



u^+Vu m -V(pdxdt = lim ( / ~u~- - Vw™ ! ■ V^dxdt ) = 0. 



This holds for all test functions <p, so the proof is complete. □ 

4. Stability of Dirichlet problems 

In this section we prove a stability result for Dirichlet boundary 
value problems in bounded domains. We use Theorem I3.1[ a local L°° 
estimate (Proposition 14.4ft . and a simple energy estimate (Lemma 14. 3ft . 
In this section, we take Q to be bounded. 

Definition 4.1. Let u G L m+1 (Q) and g G H\0,T; H\Q)). A pos- 
itive function u such that u m G L 2 (0, T; if 1 (fi))is a solution of the 
initial-boundary value problem 

{d t u - Au rn = 0, in Q T , 
u m = g, onfflx[0,T], 
u(x, 0) = u (x). 

if u m - g G L 2 (0, T; H^(Q)), and 

•/n T ^ in in 

for all smooth test functions which vanish on the lateral boundary 

of fi T . 

For the existence and uniqueness of solutions in the above sense, see 
[27| Chapter 5]. By the usual approximation argument, we may use 
test functions <p G L 2 (0,T; H^(Q)). 

Let Ui be the solution to 

( d tUi - Au^ = 0, in Or, 
(4.2) <u? = g, on^x[0,T], 

[ui(x,0) = u (x) 

in the sense of definition 14.11 We will find a subsequence of (ui) that 
converges in a suitable sense to a function u, and show that u is a 
solution of the limit problem, i.e. satisfies 

{d t u - Au m = 0, in Q T , 
u m = g, onfflx[0,T], 
u{x, 0) = Uq{x). 

More precisely, we have the following theorem. 
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Theorem 4.2. Let rrii, i = 1, 2, 3, . . be a sequence of exponents such 
that 

mi — )■ m > m c = (n — 2) + /n as i — > oo 
Let Ui be the solutions to (14. 2 p with fixed initial and boundary values g 
and uq, where 

g G H 1 ^, T; H 1 ^)), ^ G L 1+1/m ~(ft T ), and «„ 6 L m++1 (Q). 

Finally, let u be the solution to (14. 3 p itfrf/i the boundary and initial 
values g and uq, respectively. 
Then 

(1) Ui u in L q (VLt) for all 1 < q < 1 + to. 

(2) u™ 1 w m zn L s (fi T ) /or o« 1 < s < 2k, where 

, 1 1 

K = 1 + — + • 

to mn 

(3) Vu™ 1 V« m weakly in L 2 (Q T ). 

We need an energy estimate for establishing that the limit function 
attains the right boundary values in Sobolev's sense, and for verifying 
the local boundedness assumption in Theorem 13. II To derive it, we use 
the equation satisfied by the mollified solution u*\ 

(4.4) / p^+V{u m )*-Vpdxdt = [ Uo(x) (- [ (pe-'^ds) dx 
Jn T ot Jn \°~ Jo J 

This is required to hold for all test functions ip G L 2 (0, T; Hq(Q)). 
The equation (14. 4p follows from (14. ip by straightforward manipulations 
involving a change of variables and Fubini's theorem. 

Lemma 4.3. Let u be the weak solution with boundary values g and 
initial values u . Then 

(4.5) ess sup I u m+1 (x,t)dx+ [ \Vu m \ 2 dxdt 

o<t<r Jn Jn T 



dg 



Ot 



T 

1+1/m 

dxdt+ I u (x) m+1 dx 



Proof. We test the regularized equation (14.41) with p = u m — g. This 
yields 

r du* f 

/ -7r{u m - g)dxdt+ / V{u m )*- V{u m -g)dxdt 
Jn T Jn T 



«oW I i I {u m - g)e- s/a ds ) dx. 
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We rearrange this to get 
du* 



(4.6) 



at 



■u m dxdt 



V(u m Y ■ Vu m dxdt 



r du* c 

\ — #dxdt + / V{u m Y ■ Vgdxdt 

g)e- s/u ds^ dx. 



+ / u {x) I - 



For the first term on the left, we get 



du* 
~dt 



u m dxdt 



> 



— >• 



du* 
~dt 
du* 
~dt 



u(x, T) 



{u*) m dxdt + 



du* 
~dt 



Ut 



dx dt 



(u*) m dxdt 



u*(x,T) 



m+l 



m + 1 



■ dx 



m+l 



■ dx 



m + l 

as a — > 0. Here we used (12. 8p to get the inequality. For the second, it 
suffices to note that 



V(u m Y ■ Vu m dxdt 



n T 



Vw m | 2 dxdt 



as a 0. 

The limits of the left hand side terms in (14. 6 j) are positive, so we 
are free to take absolute values on the right after passing to the limit 
a — > 0. In the first term on the left, we integrate by parts before taking 
the limit: 



du* 
~dt 



g dx dt 



,dg 



u* —— dx dt — y — / u— dxdt. 



dt 



dg 



dt 



We proceed by taking absolute values, applying Young's inequality, and 
taking the supremum over t. We get 

l+l/m 

<e ess sup / u 1+m (x, t) dx + c £ 



u^- dx dt 

dt 0<t<T Jq 

The limit of the second term on the right of (14. 6 j) is 



/ 


dg 




dt 



dx dt. 



Vu m -Vg dxdt. 

Here we simply take absolute values and use Young's inequality to get 



Vu m - Vgdxdt 



<e I \Vu m \ 2 dxdt + c f |V#| 2 dxdt 
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For the third term on the right of ( 14.61) . taking the limit yields 

u™ +1 (x) dx — uo(x)g(x, 0) dx; 



The first term needs no further estimations, and the second is negative 
so we may discard it. 
We have arrived at 

(4.7) f u m+ \x,T) ^ + r ^ u ^2 dxdt 



<e / \Vu m \ 2 dxdt + e ess sup / u(x,t) 1+m dx 

0<t<T Jq 



fc / u r t\x) dx + c\ I \Vg\ 2 + 



dg 



dt 



1+1/ rn 

dx dt 



To finish the proof, replace T in the above proof by a number < r < T 
such that 

[ u(x,r) m+1 dx > i ess sup / u(x,t) m+1 dx. 
Jn 2 o<t<T Jn 

This leads to an estimate for 

ess sup / u(x, t) m+l dx + / \Vu m \ 2 dxdt 

0<t<T Jq Jq t 

in terms of the right hand side of (14.71) . The claim follows by choosing 
e sufficiently small and then absorbing the matching terms to the left 
hand side. □ 

An analysis of the above proof shows that the constant may be taken 
to depend only on m + and m~, as m varies over the interval [m~, m + ]. 
Thus, in view of (14. 5 P and the assumptions on g and Mo in Theorem 
14. 2\ we see that the sequence (Vu™*) is bounded in L 2 (Qt), and by the 
Sobolev embedding, (u™ 1 *) is bounded in L 2 {VLt). 

We use the following estimate to verify the local boundedness as- 
sumption in Theorem 13.11 See |9j Proposition B.5.1]. 

Proposition 4.4. Let u be a local weak solution, and suppose that 
B p x [t - p 2 ,t ] m Q T - Then 

1 1 Y x 

ess sup u < c 4- udxdt J +1, 

B p/2 x[t -p2/2,t } \J B p x[t a -p' 2 ,t ] J 

where 

A = n(m-l) + 2 > 0. 
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As is discussed in |9j, the constants in estimates of this type are stable 
as m decreases or increases to one, but blow up as m — > m c or m — > oo. 
Hence we may again assume that the constant is independent of i, as 
rrii varies in the interval [m~,m + ]. Note also that we have written 
this estimate over standard parabolic cylinders B p x (t — p 2 ,t ), for 
otherwise the estimate would not be stable. 

The proof of our stability result is now a straightforward matter. 

Proof of Theorem \4-S\ We combine Lemmas 12.11 and 14.31 with Proposi- 
tion 14.41 to verify the local boundedness assumption (I3.2p in Theorem 
13.11 Thus the theorem yields a subsequence and a function u such that 
itj — y u and u™ { — > u m , pointwise almost everywhere in Qt- We may 
also assume that Vw™^ — > Vu m weakly in L 2 (Qt), by Lemma [4.31 and 
the pointwise convergence. 

To pass to the limit in the equations, we need a bound on U{. We 
have 

/ u\ +m ~ dxdt <c ess sup I u 1+m '(x,t)dx 
Jn T o<t<T Jn 

<c I ess sup / u] +mi (x, t) dx + 1 
V o<t<T Jn 

This together with Lemma 1*01 implies that the sequence (ttj) is bounded 
in L l+m (Qt)- By reflexivity, we may assume that u% —>u weakly. 

By applying the weak convergences, we see that u must be a solution 
of the limit equation with the right boundary and initial values. First, 
it is clear that 

u m -ge L 2 (0,T;H^(n)), 

since L 2 (0, T; Hq(Q)) is weakly closed, by virtue of being a closed sub- 
space of L 2 (0, T; Further, we have 

-u^- + Wu m ■ Vy2 dx dt = lim f -u~ + Vm™ 1 • Vy? dx dt 
«,, dt ^°% T at 



uo(x)<p(x, 0) dx 

n 

for all smooth test functions tp vanishing on the lateral boundary and 
at t = T by the weak convergences; by uniqueness of weak solutions, 
this means that u = u. 

Convergence in measure and a bound in IP imply convergence in U 
for any r < p. For the pointwise convergent subsequences, the claims 

Ui —> u in L q (Q T ), 1 < q < 1 + m, 
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and 

u™' u m in L s (n T ), \<s<2k 

follow from this. The corresponding convergence of the original se- 
quence then follows from the fact that any subsequence that converges, 
must converge to the same limit, by the uniqueness of weak solutions. 
This also holds for weak convergence, whence we get the remaining 
claim about the weak convergence of the gradients. □ 

Remark 4.5. Generalizing Theorem 14.21 to. e.g., Neumann boundary 
conditions is a matter of proving a suitable counterpart of the energy 
estimate of Lemma 14.31 Indeed, Proposition 14.41 is a purely local esti- 
mate, independent of any boundary conditions. We leave the details to 
the interested reader. 

5. Stability of Cauchy problems 

In this section, we prove stability of Cauchy problems on the whole 
space W 1 . As for Dirichlet problems, Theorem l3.1l is the key tool. Other 
results we use are the L 1 contraction property, and the L l — L°° regu- 
larizing effect 

Definition 5.1. Let /i be a positive Borel measure on R™ such that 

fi(R n ) < oo. 

A positive function u : IR n x (0, oo) — > [0, oo) is a solution to the initial 
value problem 

u t -Au m = inM n x(0,oo), 
u(x, 0) = ji 

if the following hold. 

(1) u G L°°(0, oo; L 1 (M n )), u G L°°(R n x (r, oo)) for all r > 0, and 
u m G L l (S x (0, T)) for all compact subsets S of R n and finite 
T. 

(2) u is a local weak solution in M" x (0, oo) in the sense of Definition 

(3) For all test functions ip G C£°(lR n x [0, oo)), it holds that 
(5.1) f -u^f -u m Aipdxdt = [ (p(x,0)dn 

Solutions with the above properties exist and are unique. Indeed, by 
approximating fi with nice initial data, we can construct solutions such 
that the estimates of Theorem 15.31 below hold, and these estimates give 
the properties of u in Definition 15.11 See [61 127] for the details. For the 
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uniqueness, we note that ( 15. ip and the integrability of u m up to the 
initial time imply that 

(5.2) / u(x, t)tj(x) dx — > / rj(x)dfji as r — > 0. 

for all ?7 G C^°(M n ); thus we may appeal to the uniqueness results in 
[TT| [23] . See also j6] for an account of this uniqueness theory. 

The Barenblatt solutions ( 12. 5ft and ( 12. 6 P furnish examples of a solu- 
tion to the Cauchy problem in the sense of Definition 15.11 One usually 
normalizes B m by choosing the constant C in (12. 5p . (12. 6p so that 

B m (x, t) dx = 1 

for all t > 0. With this normalization, B m is the unique solution to 
the Cauchy problem with initial trace \i given by Dirac's delta at the 
origin, as a straightforward computation shows. 

For Cauchy problems, we have the following stability result. 

Theorem 5.2. Let mi, i = 1, 2, 3, . . be a sequence of exponents such 
that 

mi — » m > m c = (n — 2) + /n as i — > oo, 
and let m be the solutions to 



(5.3) 



d t Ui - Au™ 1 =0 in W l x (0, oo) 
Ui(x,0) = ii 



with fixed initial trace fi. Further, let u be the solution to the limit 
problem 



(5.4) 



d t u - Au m = ml n x (0, oo), 
Ui(x,0) = fi 



with the same initial trace /i. 

Then for all compact sets S in R n and all finite T , we have 

(1) — ?• u in L q (St) for all 1 < q < m + 2/n, 

(2) -> u m in L S (S T ) for all 1 < s < 1 + 2/mn, 

(3) Vu^ Vu m weakly in Lf oc (R n x (0, oo)). 

The following theorem provides the necessary estimates for our sta- 
bility result. The admissible range of the integrability exponent q in 
(15. 7p is sharp for the Barenblatt solution. This can be checked by a 
simple computation. 
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Theorem 5.3. Let u > be a solution to the Cauchy problem with 
initial data fi such that 

\\H\\ = /i(R n ) < oo. 

Then the following estimates hold. 
For all t > 0, we have 

( 5 - 5 ) IK-,*)|Ui(R") < ||a*||- 

For every t > 0, we have 

(5.6) u(x,t) < c||/x|| 2/A r n/A . 
where 

\ = n {m - 1) + 2. 

The function u m belongs to L q (St) for all compact sets S in W n and 
all finite T, for 

2 

1 < q < 1 + • 

mn 

We also have the estimate 

(5.7) f u mq dxdt < c\\fi\\^ mq - 1)+1 T-^( mq - 1)+1 . 

J S T 

Proof. The inequalities (15.51) and (15. 6p are standard estimates for the 
Cauchy problem, see [6j [261 ■ The inequality ( 15. 7p is a slight refine- 
ment of the well-known fact that u m is integrable up to the initial time 
locally in space, and follows from the first two. For the reader's conve- 
nience, we present the computation here. By applying (15.61) in the first 
inequality and (15.51) in the second, we have 

I u mq dxdt <c|| y u||^ m,2 - 1) I uix^t-^-^ dxdt 

J St J St 

<c||/x||x(^- 1 )+ 1 f r^ mq ~ l) dt. 

Jo 

We may evaluate the integral with respect to time and obtain (15. 7p if 

n 

-j(mq - 1) > -1, 

which is equivalent with 

2 

q < 1 + . □ 

mn 

The sharp constants in ( 15. 6p . as m varies, are given in [26, p. 26]. The 
situation is similar to that of Proposition 14.41 the constants are stable 
as m either increases or decreases to one, but blow up as ra 4 m c 
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m 



or m — > oo. Thus we are again free to assume that the constants are 
independent of i as rrii varies in the interval [m~,m + ]. 

Proof of Theorem \5.2[ We use (15. 6p to conclude that the sequence (ui) 
is locally bounded in IR n x (0,oo). Thus Theorem 13.11 gives us point- 
wise convergent subsequences of {uj) and (u™*), with limits u and u 
respectively. 

Convergence in measure and a bound in L p imply convergence in L 
for any r < p. From this, the claims 

Ui^yu in L q (Sr), l<q<m + 2/n, 

and 

uT 1 ->• u m in L S (S T ), 1<s<1 + — 

mn 

follow easily. We use these convergences to conclude that 



/ 

Jr 



-u^f - u m Acp dx dt 

>x[0,oo) ^ 



= lim f / ~ u i-ir ~ u^Atpdxdt ) = / v 5 (a ; ) 0)d/i. 

for any test function <p G C^°(M. n x [0, oo)). The uniqueness of solutions 
to the Cauchy problem (15. 4p now implies that u — u. The convergences 
for the original sequence then follow from the fact that all convergent 
subsequences converge to the same limit, by uniqueness. □ 

Remark 5.4. Generalizing Theorem l5.2l to Cauchy problems with grow- 
ing initial data, described in e.g. Chapters 2 and 3] or [27, Chapter 
13], offers no additional difficulties. Indeed, counterparts for all the esti- 
mates in Theorem !5.3l above are available; see for instance |27] Theorem 
13.1] for appropriate replacements of (15. 5 p and (15. 6p . An estimate sim- 
ilar to (15. 7p then follows by repeating the above computation. With 
these estimates in hand, the proof of Theorem 15 . 2 1 requires virtually no 
modifications. We leave the details to the interested reader. 



6. Stability of Dirichlet problems revisited 

In this section, we present an alternative proof of Theorem 14.21 in the 
case rrii > 1. The advantage of this argument is that we avoid using 
the local L°° estimate, Proposition 14.41 

For the reader's convenience, we give the statement of the theorem 
before proceeding with the proof, althought this is essentially the same 
as Theorem 14.21 with the additional assumption rrii > 1. 
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Theorem 6.1. Let m i; i — 1, 2, 3, . . be a sequence of exponents such 
that and 

rrii > 1 and mi — >■ m as i — >■ oo. 
Let -Uj, i = 1, 2, 3, . . be the solutions to 

{d t Ui - Aup = 0, in tt T , 
u™ = g, ondQx[0,T], 
Ui(x,0) = u 

with fixed initial and boundary values g and u , where 

g G i/ 1 (0, T; and u G L m++1 (fi). 

Finally, let u be the solution to 

{d t u - Au m = 0, in tt T , 
u m = g, ondttx[0,T}, 
u(x,0) = u . 

with the same boundary and initial values g and uq. 
Then 

(1) Ui — >■ u in L 9 (f^r) for all 1 < q < 1 + m. 

(2) m™ 1 w m m L s (fi T ) /or a// 1 < s < 2k, w/iere 

n 1 1 
k = 1 + — + • 

m mn 

(3) Vw™ 1 -> Vw m weakly in L 2 (Q T ). 

The following elementary inequality is needed in the proof. 

Lemma 6.2. For positive t, we have 

\t a - f"| < c e (l + t a+£ + t? +e )\a - p\ 

Proof. This follows by an application of the mean value theorem to the 
function x \-t t x . We have 

4-f = t x log*, 

ax 

which is estimated by 

t x logt<c £ {l + t a+e + tP +£ ) 

for x in the interval (a, (3). □ 

Stability is a consequence of the following theorem. It provides a 
quantitative estimate of the difference of two solutions. 
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Theorem 6.3. Let the exponents rrii and m, and the functions m and 
u be as in Theorem \6.1\ Then 

\\u - «i|Ui+m(fi T ) < c\m - rrii\ 1/m , 

for indices i large enough, where the constant depends on the norms of 
the boundary and initial values g and Uq appearing in Lemma \4-3[ 

Proof. We aim at using the standard inequality 

(6.3) c\a-b\ 1+m < {a m -b m )(a-b) 

in combination of an application of Olemik's test function. The function 



r)(x,t) 



0, otherwise, 



has zero boundary values in Sobolev's sense on the lateral boundary. 
We test the equations satisfied by u and and substract the results. 
This leads to 
(6.4) 

u- Ui )(u m -uT l )dxdt = - [ V{u m -u^)- [ V{u m -uf^dsdxdt 

t n 2 



n 



V(u m - uT l ) ds 



o 



dx < 0, 



1 

~ ~ 2 

where we integrated with respect to t to get the last line. Thus we have 
cf \u - Ui \ 1+m dx dt < [ {u-Ui){u m -u™)dxdt 

JQ.T JQ. T 



T 



{u-Ui){u m -u^)dxdt 
[ {u- Ui ){uf l -u?)dxdt 



< [ {u- Ui ){uf l -uf)dxdt 
Jn T 

by f[3]) and (IQ)) . 

We proceed by an application of Holder's inequality, and get 

l/(l+m) 

(u-Ui)(uf* -<)dxdt < ( / \u- Ui \ 1+m dxdt ' 
n T \Jn T 

m/(m+l) 

x ( / \uT' -uT\ {m+1)/m dxdt ' 



STABILITY OF SOLUTIONS TO NONLINEAR DIFFUSION EQUATIONS 25 

Then we apply Lemma [6.21 inside the second integral, and get 

k w< -<\<c(l + «r +e + uT +£ )\m - m;|. 
The choice of e will be made later. Thus 

m/(m+l) 



\U. 



M_ u m\(m+l)/m dx df 



T 



m/(m+l) 

<c( / (l+M™ I+£ + < +£ ) (m+1)/m dxdt) Im-rn^l. 



T 



We put all the estimates together, and end up with 

\\U — Wi||im+l(n T ) 



^ (l + <^ +e + ?C +£ ) (m+1)/m dxd^ 



l/(m+l) 



1 1/m 



< c I / (1 + + dx dt ] | m _ m . 

To finish, we need a uniform (in i) L 1 bound for the the functions 

(6 5) u mi(l+e/mi)(l+l/m) ^.(m/mi+e/m^l+l/m) 

at least for large i. We combine the Sobolev embedding (Lemma 12.11 
and the energy estimate (Lemma 14. 31) . and get that uj KiTni , i = 1, 2, . . . 
is bounded in L 1 (f2), where Ki is given by 

1 1 

Ki = l + ~ + 



n miii 

The desired conclusion follows by proving that we may choose e small 
enough, so that the exponents in (16. 5p to are less than iKitrii for large 
i. Indeed, we may choose i large enough and e small enough, so that 
1 + e/rrii and m/rrii + e/rrii are arbitrarily close to one, and so that Ki 
is arbitrarily close to k. Hence to satisfy the conditions 

mAl + — )(1 + — ) < 2/^777.3 and mA— + — )(1 + — ) < 
rrii m rrii rrii m 

for large i, it suffices to require that 

1 

1 + — < 2k. 

m 

A computation shows that this holds if 

n-2 

m > , 

n + 2' 

which is quaranteed by our assumption m > 1. Thus above arguments 
give the desired bound 

||it - Wj||ii+m(Q T ) < c\m - rrii\ 1/m 
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for i large enough. □ 

Remark 6.4. The reason why the above proof does not work well if 
m < 1 is the use of ( 16. 3p . A similar inequality is of course available for 
m < 1, but the exponents that come up in the course of the proof blow 
up as m — >■ 1. 

Proof of Theorem \6.1[ Theorem 16.31 implies that the sequence (ui) con- 
verges in L 1+m (f2r) to u, the solution of the limit problem. For a 
subsequence we get pointwise a.e. convergence, and an application of 
Lemma [2751 then yields the pointwise convergence u™* — > u m . After the 
pointwise convergences are available, the rest of the convergence claims 
are established as in the proof of Theorem 14.21 □ 

Acknowledgement. The author would like to thank professor Philippe 
Laurencot for pointing out the references [21] and |22| . 
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